Spin dynamics of possible density wave states in the pseudogap phase of the high temperature 

superconductors 
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In a recent inelastic neutron scattering experiment in the pseudogap state of the high temperature supercon- 
ductor YBa2Cu306.6 an unusual 'vertical' dispersion of the spin excitations with a large in-plane anisotropy 
was observed. In this paper we discuss in detail the spin susceptibility of the singlet d-density wave, the triplet d- 
density wave, as well as the more common spin density wave orders with hopping anisotropies. From numerical 
calculations within the framework of random phase approximation, we find nearly vertical dispersion relations 
for spin excitations with anisotropic incommensurability at low energy u> < 90 meV, which are reminiscent 
of the experiments. At very high energy uj > 165 meV, we also find energy-dependent incommensurability. 
Although there are some important difference between the three cases, unpolarized neutron measurements can- 
not discriminate between these alternate possibilities; the vertical dispersion, however, is a distinct feature of all 
three density wave states in contrast to the superconducting state, which shows an hour-glass shape dispersion. 



I. INTRODUCTION 



The pseudogap state of high temperature superconductors 
has been studied with numerous experimental tools, yet its 
origin is not resolved.^ One view proposes that the pseudo- 
gap state is a particle-hole condensate, a density wave. Of all 
such states that break translational symmetry and has strong 
momentum dependence of the type d x i_ y i, two candidate 
density wave orders that can couple to inelastic neutron scat- 
tering have been proposed: the singlet d x 2 _ y 2 -density wave 
(sDDW)^ corresponding to angular momentum I = 2 but a 
spin singlet, and the spin density wave order (SDW); in the 
general classification of density wave orders^ the latter corre- 
sponds to I = but a spin triplet. In addition to the sDDW or- 
der, its triplet counterpart 4 (tDDW) iad x 2_ y 2, where a = ±1 
corresponding to up and down spins with the z axis as the 
axis of spin quantization has also interesting properties and 
deserves more attention* 4 - Recently, Fujimoto proposed that a 
triplet d-wave particle-hole condensate may be realized in the 
hidden order state of the URu2Si2 system^,. Since high-T c su- 
perconductors have a rich phase diagram, which hosts many 
possible competing orders, it is both important and interest- 
ing to examine the properties of various density wave order 
parameters of higher angular momentum. In this paper we 
discuss the three order parameters mentioned above. In addi- 
tion, we note that a singlet chiral id x i_ y i +d 2 .j,-density wave£ 
as well as iad x 2_ y 2 +d xy density wave states with interesting 
topological properties have been explored.- Owing to limita- 
tions of space, we do not discuss these order parameters here. 

Inelastic neutron scattering can directly probe magnetic ex- 
citations. The scattering cross-section is proportional to the 
magnetic structure factor, which is proportional to the imagi- 
nary part of the dynamic spin susceptibility via the fluctuation- 
dissipation theorem 8 . Thus, a calculation of the spin suscepti- 
bility will provide a link between theoretical models and neu- 
tron scattering experiments. 

In particular, we want to address a recent experiment in 
underdoped YBa2Cu306.6- The most striking aspect of this 



experiment is a vertical dispersion relation of the spin excita- 
tions with a large in-plane anisotropy in the pseudogap state 
in contrast to the 'hour glass' dispersion observed in the su- 
perconducting stated The qualitatively different behavior be- 
tween the superconducting and the pseudogap states suggests 
different mechanisms. Motivated by the experimental obser- 
vations, we study the spin susceptibility of the three density 
wave orders mentioned above with hopping anisotropy, which 
breaks C4 rotational and mimics an 'electron nematic' stated 
In a phenomenological model, we set the hopping terms to 
be anisotropic along a- and 6-axes, and study the energy- 
momentum dispersion relations of the dynamical spin suscep- 
tibility. 

The structure of this paper is as follows: in Sec. |Il] we 
sketch the calculation of the spin susceptibility, and discuss 
the numerical results of the sDDW order. In Sec. [HI] we dis- 
cuss the numerical results of the tDDW order. In Sec. IIVI we 
also discuss the numerical results of the SDW order. To make 
the paper succinct and more accessible, the explicit forms of 
the spin susceptibility are shown in AppendixlAl 



II. SPIN SUSCEPTIBILITY: SINGLET DDW 

In this section we set up the calculation of the spin suscep- 
tibility using sDDW as an example. In the following sections 
we will give the results of the other order parameters. To cap- 
ture the in-plane anisotropic feature of the pseudogap state in 
the neutron scattering experiment, we consider the sDDW or- 
der with anisotropic hopping terms. In the momentum space, 
the order parameter can be written in terms of the fermion op- 
erators as 



(1) 



with Wk = 2^ [cos(k x a) — cos(k y b)], where a and b are lat- 
tice constants. For orthorhombic YBa2Cu306.6, a, and b are 
unequal, but the difference is very small, (a = 3.82A, b = 
3.872.) 
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The two-dimensional mean field Hamiltonian will be 

U sDDW = ^2 E ( efcC /U c fc,^ + € k+Qc\ +Q lJ C k+ Q^ 
<y k 

+iW k c\ a c k+Q , a + /i.e.) , (2) 

where the summation is over the reduced Brilloin Zone (RBZ) 
bounded by (k y b) ± (k x a) — ±7r, Q = (ir/a, n/b) is the 
nesting vector, and e k = c\ k + £2k withii 

ei fc = -2t [(1 + r) cos(k x a) + (1 - r) cos(fc^)] , (3) 
62fc = 4t'cos(fc 2: a) cos(kyb) — jj, 

-2t" [(1 + r) cos(2/c x a) + (1 - r) cos(2k y b)} . (4) 

For r 7^ 0, we have anisotropic hopping terms which 
breaks four-fold rotational symmetry. Note that although the 
anisotropy also modifies the next nearest neighbor hopping, it 
is simply a parameter and is defined as t' in our model. The 
eigenvalues of the Hamiltonian are \ k _± — ^2k ± E k with 

The one-loop spin susceptibility in the momentum and Mat- 
subara frequency space is defined as, N being the number of 
lattice sites, 

XSW,«*0 = dre^{TS\{r)Sl ql ), (5) 

where i,j = x, y, z, r is the imaginary time, T is the time- 
ordering symbol, and the spin operators are 

k,a,p 

Here a a p are the Pauli matrices with spin indices a and j3. 
We can define the longitudinal and the transverse suscepti- 
bilities as Xo z (1t1' ^ u} ) an d Xo~{Qi l' > w )> respectively, with 
= Sq±iSv and iu n — > oj + id. For unpolarized measure- 
ments, the scattering intensity, /, contains both the spin-flip 
and the non-spin-flip channels, / oc (\ zz + )/3. How- 
ever, in this paper we will present the longitudinal and trans- 
verse susceptibilities separately so that it can provide more in- 
formation about the polarized neutron scattering experiments, 
which may be achieved in the future. 

For the sDDW order, x^ilil' = 2%o (g,g',a;) be- 
cause up-spin and down-spin parts of the Hamiltonian are 
identical. The explicit form of the one-loop susceptibility is 
shown in Eq. (1AU IA4b . and we apply random phase approx- 
imation (RPA) to obtain the RPA susceptibility as shown in 
Eq. dA5HA6b in Appendix [A] For illustrative purposes, we 
set t = 0.15 eV, t' = 0.32t, t" = 0.11/ , W a = 0.65i, 
r = —0.1, and ksT = 0.05t. The chemical potential is 
set to /i = — 0.805t in order to obtain a hole doping level of 
nh ~ 10.07%, approximately the doping level in the experi- 
ment. Other similar choices of the parameters will not change 
the conclusions. 

In Fig. [Tj the constant energy cuts of the imaginary part of 
the transverse spin susceptibility along a*-axis for cj < 0.61 
are plotted. The results along &*-axis are similar and are 



not shown here. Away from Q = (ir/a,ir/b), the mag- 
netic excitations are peaked at the incommensurate positions 
(q x a, q y b) = (tt ± S a , ir) and (it, tt ± 6b), where we define the 
incommensurability 5 a and 6b along a*- and 6*-axes, respec- 
tively. From the numerical results, one finds that 6 a and 6b 
are weakly energy dependent, similar to the inelastic neutron 
scattering experiment 9 . Furthermore, a prominent anisotropy 
in the incommensurability 6b < 6 a can be seen. With the hop- 
ping anisotropy r = —0.1, we obtain 6 a w (0.30 ± 0.01)7r 
and 6 b ~ (0.235 ±0.015)tt, which gives 6 b /6 a w 0.78, which 
would be again similar to 6b/ 6 a ~ 0.6 reported in the neutron 
scattering experiments^ 

One may further adjust the parameters of this model to fit 
the experimental data, but that is not the goal of this paper. 
We have varied the chemical potential, /i, to check how the 
dispersion relations vary with hole doping; results for different 
doping levels are qualitatively similar. In the doping range 
8% < rih < 20%, there are always weakly energy-dependent 
incommensurate excitations, and the incommensurability 6 a 
and 8b increase with increasing doping level nh as shown in 
Fig. H] 




FIG. 1. (Color online) Constant energy cuts of Imx^p^(g,i^) along 
a*-axis when q y — n/b and O.lt < ui < 0.6i for the sDDW or- 
der. The weakly energy-dependent incommensurate peak positions 
are marked with red dashed lines. The results of ImxfjpA^) [J ) are 
similar and omitted. 



Note that hopping anisotropy is not necessary for the exis- 
tence of the nearly vertical dispersions. To demonstrate this, 
the numerical results with isotropic hopping are plotted in 
Fig. [3] Here r is set to 0, fj, = — 0.806t, and the hole dop- 
ing level is nh = 10.03%. All the other parameters are the 
same as in Fig.Q] One can still find nearly vertical dispersions 
with incommensurability 5 a (0.255±0.015)7r even without 
the hopping anisotropy^ 

The neutron scattering experiments show vertical disper- 
sions in the energy range 30 meV < lo < 60 meVr- and 
the numerical results exhibit a nearly vertical dispersions up 
to uj < O.Qt = 90 meV with the chosen parameters, which are 
similar to experiments. It is interesting to see how the excita- 
tion peaks evolve at higher energies, so in Fig. [4] we present 
the numerical results along the a*-axis for 0.7t < tu < lAt, 
where all the parameters are the same as in Fig. [T] except for 
the energy u>. The results along 6* -axis are again so similar 
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FIG. 4. (color online) Constant energy cuts of \mx\pA (<Z> u ) along 
FIG. 2. Doping-dependence of incommensurability S a and S b . Here a »_ axis when & = ^ / h and 7t < u < x 4t for the sDDW order . 

ii is adjusted to obtain different doping levels, and all the other pa- xhe energy-dependent incommensurate peak positions are marked 
rameters are the same as in Fig.Q] with blue dashed lines 



b 




FIG. 3. (Color online) (Color online) Constant energy cuts of 
bnXiip A (q, w) along a*-axis when q y — ir/b for the sDDW order. 
Here r = 0, fj, = — 0.806f, and all the other parameters are the same 
as in Fig.fj] 



that they are not shown here. In Fig. |H one finds that the 
high energy spin excitations are strongly energy dependent. 
The incommensurate peaks move toward q = Q in the range 
0.7t < u> < 0.9t, and eventually disappear atw « l.Ot, where 
the intensity around q = Q is enhanced. When u ps l.lt, a 
central peak emerges at the commensurate position q = Q. 
As the energy is further increased, the central peak splits into 
to two peaks deviating from Q with incommensurability S' a 
and S' b , which are marked by dashed lines. Unlike the low- 
energy incommensurability S a and <5b, 3' a an d &b are energy- 
dependent and increase with increasing energy. Note that to 
observe 8' a and 6' b , the neutron scattering experiment needs to 
be performed with very high energy (u > l.lt = 165 meV). 

The reason for the unusual vertical dispersions at low en- 
ergies and a different behavior at high energies can be under- 
stood by examining the imaginary part of Eq. iA3i . In this 
equation, the first two terms are interband contribution aris- 
ing from the scattering from the upper band (e2fc + to the 
lower band (e2k+q — Ek+ q ), and the scattering from the lower 



band (e2k — Ek) to the upper band (€2k+q + Ek+ q ). The last 
two terms, on the other hand, are intraband scattering. For the 
purpose of illustration, an example of the band structure and 
the scattering process is plotted in Fig. [5] where the interband 
and intraband scattering are shown with arrows. 

A t , ± +MeV) 




(0,0) 



FIG. 5. (Color online)Energy spectrum (\k,± + n) of the sDDW 
system as (k x a,k y b) goes along the route: (0,0) — > (7T, 0) — ¥ 
(it, 7r) — > (0, 0). The blue (red) arrows indicate the interband (in- 
traband) scattering, and the brown line is the chemical potential /i. 
The parameters are the same as in Fig. [TJ 



The interband and intraband terms of Eq. dA3l ) for O.li < 
u) < 0.6t are plotted in Fig. [6] and Fig. [7] respectively. The 
results for higher energy 0.7t < oj < lAt are not shown be- 
cause they are very similar. From Fig. [6] and Fig. [7] one finds 
that the intensity near q = Q is mainly from the contribution 
of the interband terms, whereas the contribution of the intra- 
band terms arise when q is away from Q. From Eq. dA3b . 
we can see that at q = Q, the intraband terms vanish and only 
the interband terms contribute, leading to magnetic excitations 
peaked around ui w l.lt. In the vicinity of q = Q, interband 
terms still dominate, and we may expand them to first order in 
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5q = \q — Q\ and obtain 

IT J2k ["f(£2/c ± E k ) - n F (e 2k+q T E k+q )] x 
S(uj - e 2k T E k + e 2 k+ q T E k+q ) 
- fE k [ n F^2kTE k )~n F (e 2k ±E k ) 

+^§^\E^ 2kTE Jk(e2kTE k )-Sq x 

5(u T 2E k + V fc (e 2 fe T E k ) ■ Sq), 

which will be peaked at Sq = (±2E k — a;)/ [Vfc(e 2 fc T -Efe)] • 
However, for low energies, the energy conservation condi- 
tion cannot be satisfied unless E k is very small, which di- 
minishes the difference between the Fermi functions and thus 
suppresses the intensity. Therefore, there is no enhanced peak 
in the vicinity of q = Q for low energies. For higher ener- 
gies, the energy conservation factor will be satisfied, and the 
intensity at the incommesurate positions (5' a and 5' b ) will be 
enhanced and the excitation peaks can be seen as uj > l.lt in 

Fig. a 

In contrast, away from q = Q, the intraband terms dom- 
inate. The peak positions of the energy conservation fac- 
tor, 8(uj - e 2 k T E k + e 2k+q ± E k+q ), move away from Q 
with increasing u>. On the other hand, the coherence factor 
[1 + (eifc£ifc+? + WkWk+ q ) I '{E k E k+q )\ vanishes at q = Q 
and develops with increasing \q — Q\. For the chosen pa- 
rameters, the energy dependence of these two opposite effects 
almost cancels out in the energy range < uj < 0.6t, lead- 
ing to the weakly energy-dependent positions of local max- 
ima {6 a and <5;,) as in Fig. [7] Such a dispersionless feature 
is sensitive to the parameters because it depends on whether 
the contribution of the intraband terms overcomes that of the 
interband terms away from Q. The nature of the excitation 
peaks due to the interband terms is distinct from the intraband 
terms. The dominant contribution of the interband terms are 
determined by the energy conservation factor and the Fermi 
functions, leading to sharper excitation peaks at (tt ± 5' a ,ir) 
and (tt, tt ± 5' b ), whereas the intraband terms also depend on 
the coherence factor, resulting in relatively broadened local 
maxima instead of sharp peaks at (tt ± 8 a , tt) and (tt, tt ± 5b). 
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FIG. 6. (Color online) Constant energy cuts of the interband terms of 
Imxdiag(q, u) in Eq. l lA3t along a*-axis when q y = jr/b for O.lt < 
uj < 0.6t. 




FIG. 7. (Color online) Constant energy cuts of the intraband terms of 
lmxdiag(q, u) in Eq. l |A3| l along a*-axis when q y — tt jb for O.lt < 
uj < 0.6t. 



III. THE TRIPLET d-DENSITY WAVE ORDER 

We now consider the tDDW order, and choose the spin 
quantization axis to be the z-axis without any loss of general- 
ity, that is, 

( C l+Q,a c k,p) °< ■ ^otfiWk = i{z ■ S al3 )W k . (7) 
The tDDW mean field Hamiltonian is therefore 



H 



tDDW 



- \ e ^ C k.a C k,<y + e k+Qc\ +Q <J C k+Q , a 



a k 



-i(jW k c\ a c k+ Q^ + h.c. 



(8) 



which has the same eigenvalues as the sDDW Hamiltonian. 
The explicit form of the one-loop and RPA susceptibilities are 
given in Eq. jA7MA12b in AppendixlAl 

The constant energy cuts of the imaginary part of the spin 
susceptibility of the tDDW order along a* -axis are shown in 
Fig. [8] The hopping anisotropy r is set to for simplicity 
and the parameters are the same as in Fig. [3] The longitudinal 
susceptibility behaves similar to the sDDW order whereas the 
transverse susceptibility is significantly different in the vicin- 
ity of q = Q. In comparison with the sDDW order, the in- 
tensity of Imxjt.p^q, uj) of the tDDW order is suppressed 
in the vicinity of q = Q. The intensity exhibits a V-shaped 
curve around q = Q at w = O.lt, which evolves gradually 
to a U-shaped curve at u) = 0.6t. Here we can also see the 
nearly vertical dispersion of the incommensurate spin excita- 
tions S a w (0.255 ± 0.015)7r. Notice that for unpolarized 
measurements, with / oc (\ zz + )/3, there will still be 
the vertical dispersion away from q = Q. 

The difference between the sDDW and tDDW order is that 
in XfiZ (q, of the tDDW order, Eq. ([A9i the W k W k+q term 
of the coherence factor changes sign and reduces the inter- 
band contribution. As a result, the intensity in the vicinity of 
q = Q is suppressed. The significant difference between the 
transverse and the longitudinal susceptibilities should permit 
one to distinguish the singlet and the triplet orders in spin- 
polarized measurements. On the other hand, the sign change 
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of WkWk+q does not affect the intraband terms as much as 
the interband terms, so the nearly vertical dispersions due to 
the intraband contribution can still be seen away from q = Q. 



Im *RPA(<?.o9v=-^) 

b 




lm xi? A (<i.x,qy=-M 




with those for the tDDW order in Fig. [8] In addition to this 
interchange, there is also a difference in the intensity around 
q = Q between tDDW and SDW, which could be observed 
if spin-polarized experiments with high resolution could be 
achieved, although one cannot be sure because of the non 
universal nature of this difference. Away from q = Q, we 
can also see the vertical dispersions of the incommensurate 
spin excitations with 5 a « 0.287T. Again, for unpolarized 
measurements, there will still be the vertical dispersion away 
from q = Q. 

To understand the swap of the susceptibilities between 
tDDW and SDW, we should compare Eq. (|A3t and Eq. JA9t 
for the tDDW with Eq. (lABI l and Eq. dA16l i for the SDW; 
we can see that at q = Q, W k W k + q = -W| in tDDW, and 
this leads to a minus sign, while in SDW does not. There- 
fore, the form of the coherence factors of SDW is opposite 
to tDDW in the vicinity of q = Q. As a result, the intensity 
of Imx~ftpA(q,u}) for SDW in the vicinity of q = Q is en- 
hanced due to the dominant interband contribution, whereas 
the intensity of lmxpp A (q, ui) is suppressed in the vicinity 
of q = Q. Thus, the difference in coherence factors leads 
to the "interchanging" behavior between tDDW and SDW; 
the different momentum dependence of the order parameters 
also leads to distinct momentum dependence around q = Q. 
Away from q = Q, on the other hand, both Imx^p^ (q, uj) 
and lmxp PA (q, uj) show vertical dispersion relations due to 
intraband contributions. 



FIG. 8. (Color online) Constant energy cuts of Imxlfp a(<1> uj) (up- 
per) and Imx^p A {q, (lower) for the tDDW order along a*-axis 
when q y — n/b. The parameters are the same as in Fig. [3] 



IV. THE SPIN DENSITY WAVE ORDER 

Finally, we also consider the SDW order, which has the or- 
der parameter 

( C fc+Q,a C M> K ( £ • °W0A s . (9) 

The SDW mean field Hamiltonian will be 

UsDW = ^2Y1 ( £ fc C iL C fc><7 + e k+Qc\ +Qa C k+Q}(J 
cr k 

+aA s c\ t(T Ck+Q,a + h.cA , (10) 
where the eigenvalues now become A? , = €2k ± E§ with 



b 



v fc,± ^ ^ ^k 
A?. The explicit forms of the spin susceptibil- 



ities are given in Eq. ( 1A131IAP7I ) in AppendixlAl 

The constant energy cuts of Imx^p^iq, w) and 
\mXff~ p ' A (q,ijj) for the SDW order along a*-axis are plotted 
in Fig. |f] Here we set the SDW gap to be A s = 0.65i and 
/t = -1.026*. The hole doping level is n h = 10.02%. The 
results are interesting: Imx| f 2 p j4 (g, uS) and lmxp]p A (<?, uj) 
for SDW order seem to be 'interchanged' in comparison 




FIG. 9. Constant energy cuts of lmXiipA(<li ljJ ) (upper) and 
lmx^ip A {q, uj) (lower) for the SDW order along a*-axis when q y = 
n/b. Here A s = 0.65t, /i = — 1.026f, and the other parameters are 
the same as in Fig. [8] 
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V. CONCLUSION 



Appendix A: The explicit forms of spin susceptibility 



In conclusion, we have attempted to provide an explanation 
of a recent neutron scattering measurement in an underdoped 
high temperature superconductor, which point to the fact that 
the pseudogap state is not a continuation of the superconduct- 
ing state below T c . The salient feature is a vertical dispersion 
seen above T c , in the spin excitations as opposed to an hour- 
glass shape dispersion seen below T c . 



Although couched in the language Hartree-Fock theory 
augmented by RPA, a thorough analysis of the properties of 
various alternate order parameters should be a useful guide. 
We also checked a band structure to contain electron pock- 
ets as well, but the robust aspects of the conclusions were 
unchanged. The vertical dispersion feature appears to per- 
sist in the doping range 8% < n/, < 20%. At higher en- 
ergies, we find energy dependent incommensurability due to 
the interband contributions. We also contrast the spin dy- 
namics of the tDDW and SDW orders, which exhibit differ- 
ent features around q = Q, which could in principle allow 
one to identify the spin nature of the underlying phase in a 
spin-polarized neutron scattering experiment with high reso- 
lution. The transverse and the longitudinal spin dynamics are 
interchanged between SDW and tDDW. In principle, a whole 
class of higher angular momentum particle-hole condensates 
are possible. Experimental evidence of these order parameters 
should be a major step forward. The tDDW is such an uncon- 
ventional hidden order that its discovery would be of great 
importance. Note that tDDW is even invariant under time re- 
versal. 



ACKNOWLEDGMENTS 



This work is supported by NSF under Grant No. DMR- 
1004520. 



In the density wave systems we considered above, the 
Green's functions form matrices 



G a (k,iLo m ) 



g a (k,k,iui m ) g a (k,k + Q,iu) m ) 
g a (k + Q,k,iLo m ) g a (k + Q, k + Q,iu m ) J' 



where 



g rT (k,k',iuj m ) 



dre^(Tc k Ar)clj 



The one-loop spin susceptibility also has diagonal and off- 
diagonal terms 

Xo'fo?')"^) = Kq'Xti%(q,iUm) + <W+QXrfT 0?, *w m ), 

where the subscripts 'diag' and 'off refer to the diagonal and 
off-diagonal terms of the one-loop spin susceptibility, respec- 
tively. 

With a quadratic Hamiltonian, these terms can be written 
in terms of the Green's function matrices by applying Wick's 
theorem, and we have 

xlSgfo *w m ) = — ^2 Tr[G a (k + q,ie n + iu m )G^(k,ie n )], 

k,n,a 



k,n,a j^l 



Xdiag (?. *w m ) = -^^2Tr[G r {k + q,ie n + iui m )Gi{k,ie„)], 

" k,n 

xta ~{q,iu m ) = jjy^2^2[Gt(k + q,ie n + iu} m )Gi(k,ie„)] 



where TV is the trace, and G a (k,ie n ) can be obtained from 
the Hamiltonian. 

For sDDW, the up-spin and down-spin components are 
identical. For a =t or \,, we have 



G a (k,ie) 



i( + fit - £2fe iW k 



(if- - i-2k) 2 - El \ -iWk ie - eifc - £2fc 
Therefore, we have 



Xo 2 (<7,<7',w) = <*<,,<,' Xdiagfow) + <W+QXaff(g,w), (Al) 



Xo fog')") = 7rX(Tfo</^), 



where 



(A2) 



Xdiagfow) = 



*£ fc (i- 
-*£*(n 



£lk£lk + q + W k W k + q s 

E k E k+q i 

. tlk£lk + q+W k W k + q 
E k E k + q 



nF{e2 k +E k )-n F (e 2 k + q-Ek + q ) , n F (e 2k -E k )-n F (e2 k +q+E k + q ) 



-(■2k — E k +€ 2 



-Et 



-2k + q~ 

nF{e 2k +E k )-n F {<i2k + q+E k + q ) 
u-e2k-E k +62k + q +E k+ q+iS 



U-t2 k +E k +t2 



a+E k+ „+tS 



7-F(£2 k —E k )—n F (E2k + q—E k + q ) 

u-e 2k +E k +e2 k +q-E k + q +iS 



(A3) 
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E k E 



k t q 



n F {t2k+E k )-n F (<i2k + q-E k+q ) _ n F {t 2k -E k )-n F (t 2k+ q + E k+q ) 



'-tzk-Ek+tzk + q-Ek+q+ii 



u-e 2k +E k +t 2k + q +E k + q +i8 



■n F (e 2k +E k )-n F (e 2k + q +E k+q ) 
u — £2k —E k +e 2 k+q+E k+ q+iS 



n F {<i 2k ~E k )-n F (e 2k + q -E k+q ) 
u-£2k+E k +e 2k + q-E k + q +i5 



(A4) 



where np(E) is Fermi-Dirac distribution function, and 6 is 
set to 0.06< for the numerical calculation in order to obtain 
smooth curves. Applying random phase approximation, we 
obtain the RPA susceptibility^ 



Uxl z {q\,q' » 

-+- i ' \ Xo"~(?,9i> w ) 



(A5) 



(A6) 



where Xo z (?, 9', w ) anc ^ Xo~~ (?, 9', w ) ^ me two by two ma- 
trices from Eq. (lAlb and Eq. (IA2t . respectively. For the nu- 
merical calculation, we set U — Wo — 0.65t and compute the 
imaginary part of the diagonal terms of the RPA susceptibility 



(q = q'l 

For the tDDW order, the Green's function matrices become 



G a {k,ie) 



«e + eife - £2fc 



icrWk 



(ie - €2k) 2 - E? \ -ivWk it - e lk - e 2k 



where a = +1 for up-spin and a = — 1 for down-spin, and 
the spin susceptibility will become 

x z z (q,q',^) = <Wx£g(9iw)i ( A? ) 

Xo~(q,q',u) =6 qiq >Xfa e (q,u) + 6 q , q > + QXtfT(q,u), (A8) 
where x$L(<l> w ) i s me same as Xdia g (<Z, w) in Eq. ( |A3t , and 



X<L(<?^) 



1 v-» (-1 _ e lfc e lfc+g -W fc W" fc + q ■, 

ei k e lk + q -W k W k+q ■ 
E k E k+a 



n F (t 2k +E k )-n F (t 2k + q -E k+q ) ^_ n F (e 2k -E k )-n F (e 2k+q + E k+q ) 



2N Sfe(l 



"F(e2fc+-Efc)-"F(e2fc+ g +-Efc+ g ) 
u-e 2k -E k +e 2k+q +E k + q +i5 



J — £ 2k J rE k +e 2k+q +E k -t- q +i5 

•n F (e 2k —E k )-n F (e2 k + q -E k+q ) 
u — £ 2k +E k +c 2k+q — E k + q +i8 



(A9) 



— i f e ikW k j r g+ei k ^ q W k 

-2N l^k\ E k E k + q 



n F {£2 k +E k )—n F (e 2k + q —E k + q ) 
u-e 2k -E k +e 2k+q -E k + q +i5 



n F (t2 k — E k ) — n F (e 2k + q +E k+q ) 



n F (e 2k +E k )-n F (e 2k+q +E k + q ) 
w—£ 2 k—E k +e 2k + q +E k + q +i8 



£ 2 k+E k +e 2k + q +E k+q +iS 

n F (t 2k -E k )-n F (t 2k+q ~E k + q ) 
u — e 2k +E k -\-€2k+ q —E k + q -\-id 



(A10) 



r 



The RPA susceptibility of the tDDW order will be 



xSm(<?,?', w ) = 



Xo z (q, q',^) 



l-U X z z (q,q',u) 
^ V J -^-(«i J 3 , J w)' 



(All) 
(A12) 



where Xg z (g, q', w) is from Eq. (IA7b and %q (q 7 q' 7 uj) is a 
two by two matrix from Eq. 



For the SDW order, the Green's function matrices become 

G a {k,ie) = 



1 



crA s ie - £i fc - e 2/c 



( l£ - £2fc ) 2 -(£;f) 2 
The longitudinal and transverse spin susceptibility are 

X z z {q,q',uj) = <WXd?ag(3,w), (A13) 

Xo~0?,?',^) = Kq'X&agiq,") + <W+QxJ fi r (<7,^), 

(A14) 

where x|L (?, w), xj[" (g, w), and x+ ff ~ (g, w) now become 
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Xdfag(<7>^) 



Xdia a (<?> w ) 



e lfc e lfc+g 



n F (e 2k +El)-n F ^ 2k+q -El +q ) 
ui-e 2k -E^+e 2k+q -E^ + +iS 



n F (e 2k -E k )-n F (e 2k+q + E k+q ) 
u-e 2k + E% +e 2k+q + E* +i s 



1 l-\ I EltElt+l+^il 



"F(e 2 fc+-Ef )-n F (e 2k + q +E^ +q ) n F (e 2k -E^)~n F (e 2k+q -El +q ) 



2N J2k(^ 

fe + g 



"F(e 2 fc+gf )-n F (e 2 fc+g-gfc +g ) _ 
u-e 2k -E^+e 2k+q -E^ +q +iS 

n F (e 2k +E^)-n F (e 2k+q +E^ +q ) 
uj-e 2k -E^+e 2k+q +E^ +q +iS 



n F (t 2k -Ej)-n F {e 2k+q +Ej +q ) 



u>-e 2k +E^+e 2k+q +E^ +q +ii 
n F (e 2k -E^)-n F (e 2k+q -E^ +q ) 



-£ 2 k+E£+e 2 



Ef E ;- 



u+e 2k +Eg-e 2k+q -E^ + +i6 



"+£ 2fc --Ef-e 2fc + ,+£f + „+i<5 



, / ^f +-Ef+, N "F(e 2fc +gg)-"F(e 2fc+ g 7 -Eg +g ) 
1 ^^+, j "+e 2fc +B fc i - £2fc+ ,+B^ + ,+ l 5 



- g fc+^+g ^ "F(e 2fc --Eg)-»F(e 2fc + g+gfc + g) 

~Ej7~+iS~ 



EfEf +q ' u,+e 2k 



V 

1 "+*—E k *-£2k+ q -E k i +q - 



(A15) 



(A16) 



(A17) 



The RPA susceptibility of the SDW order is in the same form as the tDDW order in Eq. (IA1 lb and Eq. (IA121 . 
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